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Abstract 

A spherically symmetric comoving fluid solution of Einstein's equa- 
tions is adapted for cosmological application by extending the geometry 
of standard FRW cosmology using a generalised curvature term. The re- 
sulting model retains many of the known cosmological properties includ- 
ing homogeneity of energy density, its relationship with internal pressure 
including equations of state, although in each case they have a gener- 
alised structure. It is shown that the adapted model does not require the 
inclusion of the arbitrary cosmological constant and the vacuum energy 
solution is discussed in its absence. The Hubble constant and deceleration 
parameter are also shown to have a form which characterises the modified 
geometry of the new model. These forms are calculated using current ob- 
servational data and show how the standard cosmological geometry can 
be amended in a way which is consistent with an observed flat curvature 
and a decelerating universe. Finally the solution is also considered in 
the context of gravitational collapse where it is shown how fluids spheres 
obeying a central equation of state can be matched to empty spacetime. 

Short title: Relativistic fluid spheres applied to cosmology and grav- 
itational collapse 
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1 Introduction 

Spherically symmetric fluid solutions of Einstein's equations have been discussed 
in depth by many authors for many reasons which include their important wide 
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ranging applications to cosmology and astrophysics. General solutions have 
often been obtained using a metric in isotropic form together with a comov- 
ing system of coordinates. From the point of view this paper the approach 
originally described by Kustaanheimo and Qvist [T], has been invaluable al- 
though although there was an earlier particular solution presented by McVittie 
. Subsequently there have other investigations which often involved Lie sym- 
metry approaches for example Stephani [3], Stephani, Wolf [4| and Wiltshire 
[S] . In addition have been many comprehensive studies of their mathematical, 
physical properties and their various interrelationships by for example Nariai [B] 
,Chakravarty [7], McVittie §1^, Srivastava [TO], Sussman [HinillB], McVittie 
[5] , Knutsen [T3] . Furthermore the host of solutions have been summarised by 
Krasinski [15] and also Stephani et al [16]. 

Nonetheless there remain many outstanding problems associated with the 
physical applicability of these solutions. For example in cosmology the simplest 
spherically symmetric solution gives rise to the Friedmann, Robertson, Walker 
(FRW) standard cosmological models which have recently required significant 
adaptation to describe emerging observational data. The evidence of very 
early epoch, inflation, the notion of vacuum energy and a negative deceleration 
parameter in a more or less flat universe (based upon type la supernova surveys 
Perlmutter, Riess, Schmidt, Garnavich & coworkers [TT] (TS] [TH [H]) have 
created many challenges for the standard model and general relativity. This 
has required the reintroduction of the cosmological constant and the resulting 
concordance or ACDM model to reconcile observation with theory by means 
of large proportions of dark energy f^A — 0.7 and cold dark matter flM — 0.23. 
Recent developments Turner & Reiss [22], Virey [25 show that even this may 
not be sufficient, in fact a bimodal model may be required to describe past 
deceleration as well as recent acceleration of the universe. In the following it 
will be seen that there is a generalisation of FRW cosmology which overcomes 
many of these difficulties without the introduction of a cosmological constant. 

Moreover this solution will be shown to simplify many difficulties in the 
description of problems of gravitational collapse where in practice it is extremely 
difficult in practice to match a ffuid sphere to Schwarzschild empty space-time 
even though the theory is well known, Bonnor & Vickers [H] . This is especially 
true in cases where it is necessary to incorporate equation of state between 
pressure and energy density p = p{p) for example at the centre of the sphere. For 
this reason the Oppenheimer Snyder model [25] essentially a zero pressure FRW 
solution has often been employed in matching problems as has been the case in 
discussions of a potential source for gravitational waves Babak & Glampedakis 
[26] or a matching problem in first order rotation Kegeles [27] . 

This paper is organised in the following way. Einstein's equations in the 
context of general spherically symmetric fluid spheres are introduced in section 
1 with a focus on the particular solution to be discussed. The basic properties 
of the solution are presented in section 3 in relation to the FRW cosmological 
models and the cosmological constant. The solution is then discussed in the 
context of vacuum energy in section 4 whilst section 5 focusses on equations of 
state and problems of gravitational collapse. Finally in section 6 the solution 
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is discussed in terms of the Hubble constant, deceleration parameter and some 
current data. 

2 Fluid spheres, isotropy, preliminary equations 

It is the intention here to consider the isotropic coordinate system for which 

ds^ = e^^dt^ ~ e^f" {dr"^ +r^dVP) dn^ ^ d0'^ + sin^ (6) dc/)^ (1) 

where A — A(r, t), fi — fi{r,t). In addition Einstein's field equations will be 
taken to be 

Gl = -SnTl n = {p + p) u^uk - dip (2) 

where is the Einstein tensor and is the energy momentum tensor of the 
fluid sphere. A comoving coordinate system will assumed so that the compo- 
nents of the velocity four vector satisfy u^Ui — 1 with = and = = u^. 
The pressure p, energy density p and mass function m, see for example, Misner 
and Sharp [2S], Cahill and McVittie [15] may be calculated using: 

8Trp = Gl 8TTp=-Gi m=^{l + e2(^-^VV?-(l + ?-A*r)^} (3) 

where the suffix r, t indicates a partial derivative. In this system Einstein's 
field equations satisfy the isotropy condition in the form 

Gt = Q = Gl Gl-Gl=0 (4) 

and so 

Prt - KfJ-t = (5) 

and ^ 

Mrr ^~ A^^ -|- A? — pZ — 2Aj./i^ — =^ (6) 

r 

Using the approach of Kustaanheimo and Qvist [1] the solutions of ^ and ([6]) 
may be expressed in terms of the function L = L{x, t) where 

1 = 6-" x = r'^ e^ = A {t) pt = -A (t) — (7) 

In this notation ([5]) dtiid ([6]) are simultaneously satisfied by 

L,, = L'F {x) (8) 

In general this equation has been studied for the mathematically tractable 
cases F {x) — {ax^ + 6x + c)~^/^ often leading to solutions involving elliptic 
functions, see for example Stephani et al |16j whilst the much simpler form 

L{x,t) = f{t)x + g{t) F(x) = (9) 
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leads to the FRW cosmological models and also the Oppenheimer- Synder model 
for gravitational collapse [25] . However the FRW solution 



L(.,t) = ^(l + ^) F(.).0 (10) 

is only a subcase of the more general class given by © which will be the primary 
subject of study below. It will be shown how this generalised class leads to a 
new class of fluid spheres possessing uniform energy density with the prospect of 
extending the class of FRW cosmologies and also fluid sources for gravitational 
collapse. 



3 The solution and basic properties 

In the context here it helpful to defines the function S (r, t) for which 

S = 



rR 



1 



aR^r"^ 



(11) 



and a solution of Einstein's equations having the general structure of ^ in the 
form: 



e-^' = - = 
S 



1 



R 



so that by © 



St _j 



R 



r^R {2aRt + atR) 



V 



4 1 



aR'^r^ 



(12) 



(13) 



where R — R{t), a = a (t) and a — a{t) and satisfy ([5]) and ^ by direct 
substitution. The Ricci scalar curvature for the 3 dimensional space in ([T]) with 
has the value 6(7 so that the solution is closed when ct > 0, flat when cr = 
or open when a < 0. 

Also note that the FRW solution is obtained as a subcase of ^ with p2|) 
and ()13|) by writing 



k_ 

i?2 



^2 



3. 

i?2 



The energy density p ^ p{t) from ([3]) is: 

87rp = 3 (a + cr) 



(14) 



(15) 



It is a purely time dependent equation and is a generalisation of the normal 
Friedmann equation when the condition (|14p is also to be included. Equation 
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demonstrates that the energy density is a composite of two functions a {t) 
and a (t) neither of which have any a priori dependence of the scale factor R (t). 
However the pressure p = p{r, t) calculated using ^ is in general a function of 
both r and R (t) and is given through 



Stt {p + p) 



{at + at) R {4 + ar^ R^) 
{r^aR^Rt - 4:Rt + r^atR^) 



(16) 



Equations ([TC]) and represent the generalised form of the two Friedmann 
equations which form the basis of standard cosmology. Subsequent analysis in 
this paper will be based upon these generalised forms. It can be shown that 
(|16p may be cast into a more familiar form to cosmologists by combining it with 
S (t) defined in (fTTj) and the derivative of (fTS]) to give 



P- 



dt 



d {pS^\ 



S{p + p)St 



Notice that when either r 
form 



dt =^ = s — 

= or alternatively <t = k jR? this reduces to standard 



Pt 



3 (P + p) Rt 
R 

In the addition the mass function has the form 

^3 



m 



^-KpS-^ inp (Rr) 



1 + 



crR'r' 



(18) 



(19) 



also part of standard cosmology and which also be used in cases of gravitational 
collapse to define a fluid sphere boundary. 

As can be seen from equation it is not possible in general to choose a 
function a {t) such that = 1 for all values of r and so define a comprehensive 
proper time variable, that is with the exception of the FRW cases given by dH]). 
However it is possible for an observer at r = k to define a local proper time 
variable by setting 



S (k, t) = K (t) 



S{0,t) = R{t) 



(20) 



and using ([T 



ait) 



Kl 



I \ 

Rt K^R{2aRt+atR) 



V 



4 1 



(21) 



In this way becomes 



(22) 



whilst the pressure from (fT7| now satisfies 



5 



p- 



d (pK^) 



Pt = 



3{p + p)Kt 



(23) 



dt dt K 

Without loss of generality in the following it will be assumed that local 
propertime is defined at the coordinate centre r = so that: 



whilst from (|23 



Snp = 3 



2Rt 



i?2 



<JtR 



3a 



(24) 



(25) 



R i?2 Rt 

Clearly in the particular case when a — k/R^ then ((24|) and ([25|) give rise 
to the Friedmann equations. Note also that on writing 



X 



(26) 



where a — a (t) and x is a constant then ([24|) and ([25]) become respectively: 

2Rtt Rt ^tR 



8np = 3[^+a]+X 



R 



_L„^_3a-x (27) 



Thus the choice of a at (|26p results in the systematic inclusion of a constant, 
X equivalent to the cosmological constant A (= — x) normally associated with 
Einstein's equations. However, in terms of cosmological applications it is better 
to consider equations (j24p and (j25p in their full generality as cr is a priori an 
undefined function of t and so has the capacity to act as a 'variable cosmo- 
logical constant'. For example, Perivolaropoulos I 30j describes the cosmological 
constant problem whereby the cosmological constant is considered to have had 
relatively large value during the early period of rapid inflation much larger than 
during the current epoch. 

Finally note that the metric ([T]) with (I12p and (|13p remains invariant under 
the transformation 

r=- R=^ (28) 

where R — R{t). In particular the transformation gives to 



ds^ = e^'^dt^ - e^f" {df^ + f^dn'^) 



where 



R 



( \ 

R^ _ r'^R{2aRt + atR) 



V 



4 1 



cjR^r^ 



a 2 



(29) 



(30) 



Thus the properties of the energy density and pressure found at (r, t) are 
replicated at (f, i) using ((28)) . In particular the properties found at r = are 
replicated at infinity. 
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4 Vacuum energy and equation of state p = 

Consider now the case in (l26l) when 



so that equation ([?/)) reduces to the vacuum energy equation of state (see for 
example, Peacock [3T] ) 

8ttp = x P=-P (32) 

The full solution is now given by equations (fT2|) and (fT3|) 



(33) 

If proper time is employed then on writing — x/3 with (|3ip then 

R 

and also 

cr = m = ^ (35) 

When X — ^ then from a Minkowski spacetime is obtained represented by 



Thus using proper time 

with ^ 

cr = 2 = (38) 

(1 + at) 

Equation (|37p is of a form described by Milne [32] and Peacock [31] . 



5 Equations of state 
5.1 Cosmological application 

Consider the cosmological case so that there is no fluid boundary and where 
it will be assumed that a fluid equation of state p — p (p) exists at r = 0. In 
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addition local proper time is used so that so that from equation (PT|) 

« = || (39) 

With K = R, equation ((23| may be solved for particular p = p{p) and the 
results which have a familiar form are summarised in Table 1. 



Equation of state 

( at r = 0) 


Pressure p 


Energy density p 

(a is constant) 


Adiabatic 

(7^1) 


p = Np^ 


i 


Linear general 
(all values of n) 


p={n- l)p 


a 


Dust 


p = 


a 


Radiation domination 
[n = 4/3) 


P 

^=3 


a 


Vacuum energy 
[n = 0) 


P= -P 


p = a 



Table 1 



Thus the function a (t) can then be determined from equation (|24p so that 

(40) 



Snp Rl 



3 i?2 

This equation determines the function a {t). In FRW cosmology a = k/R^ then 
equation is a differential equation that is solved to determine R = R {t). 

5.2 Gravitational collapse application 

In cases of gravitational collapse suppose that a fluid sphere has an equation 
of state p = p{p) defined at the centre r = which from equation pT|) gives 
results which are identical to those described in Table 1. In this context the 
adiabatic equation of state is particularly relevant to polytropic stars as has 
been described by Weinberg [53]. However it is not assumed that local proper 
time is used at r = and so equation ((2T|l does not hold. However it is 
supposed that the fluid body has a well defined boundary r = b which matches 
the Scharzschild vacuum solution and where the boundary value of the pressure 
must be zero. 

For zero pressure at r = 5 it is well known, see for example, Bonnor & Vickers 
[23] or Cahill and McVittie [IS] that the mass function (fTOj) is a constant m = M. 
Thus 

M.l^ = l!^ (41) 
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where 

Rb 



4 

so that 



i^ = 5W = ^^ (42) 
1 + 



|(.B')^0 ^ B.-lj (43) 



Equations (|^^ together with may used to show that 

4 



whilst from equation ([Tl 



K + Rbp^/^ ) (44) 



3 R'^b^K 
It follows for dni) with (gSl) that 

^ dt dt ^ ' 

which gives the expression for pressure such that p {b, t) = as required. 

Finally note that each of these fluid spheres when endowed with first order 
rotation in terms of a rotation parameter may also be matched to empty space- 
time using a general result given by Wiltshire [34] . 

6 Hubble constant and deceleration parameter 

A further understanding of the geometry of this solution and its relation to 
current cosmological observational data may be obtained by determining the 
Hubble constant and deceleration parameter for the solution. 

The null geodesic equation for an inward travelling photon as described by 
an observer at (r, t) — (0,io) using proper time is 



dt 

With the notation that 



-e^-f" with = 1 (47) 



t = tQ-T (48) 

where r is the travel time of the photon and using and equation II47P 
expressed in terms of travel time r becomes: 



dr 1 + r^aR"^ )^ _ 



dr R 



( \ 

Rt_ r^R {2aRt + atR) 



a-2 = 1 (49) 
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The equations (|^^ are solved upto and including second order terms in r by 
writing 

R = Rq- Rot + —T o- = (To - o-qt + —T 

a = ao- aoT + '^t^ (50) 

where for example R means Rt evaluated at i = to- In this way the solution the 
first of (I49p and for example the redshift z can be calculated in the usual way 

r (r) i?o = r + and ^ = Ht + H'^^^t^ (51) 

where the Hubble constant and deceleration parameter are 

rr Ro RqRo /t.o\ 

The condition that = 1 is given by 



-^0 ■ 2i?0-Ro „-Ro /roN 

ao «o = —S2 2^ (53) 



and is essentially the statement that a = {Rt/ R)^ upto and including first order 
terms in r. 

With this notation the first of with and gives 

^ ^ ~ ~3 ^ 



whilst the second of ([53]) with ([52|) together with (|T8|) results in 
i?o^o o'o 47r (po + Po) 



1 (55) 



So when p = [n — l)p then equations ((54)) and ([55]) may be taken to give 

I : — I 



Thus in the case when a — {) for all t then 

<l-'i-l (57) 

as expected. In addition for the case of vacuum energy n — Q and a is constant 
then q = —1 again as required. 
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In the light of recent observations for example WMAP [3S] evidence from 
the cosmic microwave background suggests that the universe is essentially flat 
so that in the context of this model (Tq = so that 



H'=^-^ (58) 

as is normally calculated. In addition however observations of distant type la 
supernovae suggest that q <0 and that the universe is decelerating rapidly. In 
this model (|56p can be rearranged for ctq incorporating — so that: 

^0= (^)'^'(2<z-3n + 2) (59) 

With equation ([55)1 this is also 

(To = {2q -in + 2) (60) 

Thus for example a matter dominated universe, n — 1 with q = —1/2 gives 

cro = -2H^ (61) 

Hence it has be shown how the observed data can be used to determine the 
geometric parameters [uoj o'o] which characterise the difference between the ge- 
ometry of the current model with that of standard FRW cosmology. 



7 Conclusion 

In this paper the primary focus has been on consideration of an extended ver- 
sion of the FRW solution which can be characterised by a Ricci curvature scalar 
having a value Qa{t) which is independent of the scale factor R{t). Moreover 
the purely time dependent energy density consists of the sum of two terms which 
from a geometrically point of view are the curvature and also a function a {t) 
that defines the nature of the time coordinate, for example proper time. This 
is an extended form of the Fricdmann equation and contrasts with the FRW 
solution for which the curvature and also the energy density depend closely on 
the scale factor. It has further been shown how expression for internal pressure 
now a function of (r, t) gives rise to an equation which is closely analogous to 
the second Friedmann equation which expresses the acceleration of the scale 
factor in the standard model. The generalised equations are then shown by a 
translation of the curvature function to contain naturally a constant that can be 
interpreted as a cosmological constant and can also be used to describe vacuum 
energy. However its introduction is unnecessary in the new model as it is better 
to to consider the curvature in its full generality for the purposes of cosmological 
application. For these purposes it is shown how to define local proper time in 
such a way that observer essentially see a homogeneous cosmology but with the 
inclusion of the modified curvature term. The Hubble constant and decelera- 
tion parameter are calculated in the usual way but now include a component 
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that reflects the gcnerahscd form of the curvature term. Hence when the obser- 
vational data (an essentially flat universe coupled with a negative deceleration 
parameter) is introduced the Hubble constant and deceleration parameter are 
shown to provide information about the generalised nature of the curvature and 
the geometry of the modified cosmological model. In terms of this model the 
data is not interpreted in terms of dark energy or cold dark matter. 

Finally the solution is discussed using a equation of state for an observer 
employing local proper time. Moreover the resulting solutions are also described 
in terms of problem of gravitational collapse. In particular it is shown how 
fluids including polytropes, with a central equation of state may be matched to 
Schwarzschild empty space-time. 
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